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ABSTRACT 

The problem of nonlinear development of Gortler vortices and 
interaction with TS waves is considered within the framework of 
incompressible Navier-Stokes equations which are solved by a Fourier- 
Chebyshev spectral method. It is shown that two-dimensional waves can be 
excited in the flow modulated by Gortler vortices. Due to nonlinear effects, 
this interaction further leads to the development of oblique waves with 
spanwise wavelength equal to the Gortler vortex wavelength. Interaction is 
also considered of oblique waves with spanwise wavelength twice that of 
Gortler vortices. 
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1. INTRODUCTION 

The subject of laminar/turbulent transition is of fundamental and 
practical importance in fluid mechanics. An indepth knowledge of the 
transition mechanism is needed not only for boundary/shear layer control 
(e.g., delay of transition or mixing enhancement), but also for 
understanding of turbulence. Depending on the state of the boundary layer, 
various instability mechanisms such as Tollmien-Schlichting (TS), 
crossflow and Gortler vortices may be operative. What instability modes are 
actually excited in a boundary-layer depends a great deal upon the 
particular "forcing" present. This is the classical "receptivity" problem. 
The presence of a finite-amplitude disturbance in a boundary layer could 
also lead to the excitation of disturbance modes which may otherwise be 
damped or weakly unstable according to linear stability theory. Examples of 
such interactions are the possible excitation of TS instability in the presence 
of Gortler vortices on a concave wall, fundamental or subharmonic 
secondary instability in a flat plate boundary layer or the excitation of TS 
instability in the presence of crossflow vortices on a swept wing. 

Consider a boundary-layer flow on a concavely curved plate. Here, 
counter-rotating steady Gortler vortices form due to centrifugal instability. 
If the Reynolds number is high enough, TS instability may also be present. 
Nonlinear development of Gortler vortices and Gortler/TS interaction is a 
problem of both fundamental and practical importance. The interaction 
may take place in three different situations: 

(i) when Gortler vortices are finite-amplitude and TS waves are 
infinitesimally small. 

(ii) when TS waves are finite-amplitude and Gortler vortices are 
weak. 

(iii) when both Gortler and TS waves are finite-amplitude so that 
they interact nonlinearly. 

Nayfeh (1981), using multiple- scale analysis for two-dimensional 
boundary layers, showed that a finite-amplitude Gortler vortex could interact 
with two oblique TS waves of spanwise wavelength twice that of the Gortler 
vortex. The resulting growth rates were found to be many times larger than 
the linear theory growth rate of the TS wave without the Gortler vortex. 
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Furthermore, these exponential growth rates were predicted to increase with 
increasing amplitude of the vortex. Since the amplitude of the vortex 
increases exponentially, there is essentially a 'double exponential growth' of 
the disturbances in the boundary layer. Similar results were also obtained by 
Srivastava and Dallmann (1987). The concept of 'double exponential growth' 
has also been suggested by Herbert and Morkovin (1979) for the interaction of 
finite-amplitude TS waves with streamwise vortices and by Floryan and Saric 
(1980) for the streamwise vortices interacting with Gortler vortices. Recently 
Hall and Smith (1988) have also studied Gortler/TS interaction. Such 
interactions are crucially dependent on amplitude of the primary mode. The 
present study is an attempt to determine the amplitudes of the primary 
instability for these interactions to take place. Other types of interactions 
including those which involve oblique TS waves with spanwise wavelength 
equal to the Gortler vortex wavelength are also discussed. 

In this work, the wave-interaction problem in two-dimensional 
boundary layers on a concave surface is studied numerically. The three- 
dimensional incompressible Navier-Stokes equations are solved by a Fourier- 
Chebyshev spectral method with two periodic and one non-periodic direction. 
Even though the boundary layer grows spatially, the works of Wray and 
Hussaini (1984), Zang and Hussaini (1985) and Spalart (1984) have shown 
that temporal solution of the Navier-Stokes equations with periodic boundary 
conditions in the streamwise direction is able to capture the qualitative 
features of the transition process on a flat plate quite adequately. It is 
assumed that such will be the case in the boundary layer that is considered 
here. Though the method is applicable to all the regimes of interaction listed 
above, results are presented only for the case where Gortler vortices are 
finite-amplitude but TS waves have relatively small amplitude. 

2. GOVERNING EQUATIONS AND NUMERICAL APPROACH 

We consider the problem of laminar/turbulent transition in two- and 
three-dimensional boundary layers. The governing Navier-Stokes 
equations for an incompressible fluid are: 


U t +U. VU = -Vp + vV U + F 
V.U=0 
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(2.1) 

(2.2) 


where U = (u, v, w) is the velocity vector, p the pressure, and v the kinematic 
viscosity. The term F is a forcing term. No-slip boundary conditions at the 
solid wall (z = 0) are imposed, i.e., 


U(x, y, 0, t) = 0. 

In the free stream, it is reasonable to require that 

U as z -» °o. 


(2.3) 


(2.4) 


Equations (2. 1-2.4) are solved by a Fourier-Chebyshev spectral method 
similar to that described by Malik, Zang, and Hussaini (1985). Periodicity is 
assumed in the x and y directions. The dependent variables have Fourier- 
Chebyshev series of the form 

K, K„ 

1 -^-1 

. . 2 -r-i 2 -r-> 2ni(k,x/L x ) 2ni(k y y/L y ) 

u(x, y , z, t ) = y y £ u(t)e e T n (r|) 

ky = -V=° 

2 2 (2.5) 

where Lx and Ly are the periodicity lengths in the x and y directions, 
respectively, and T n is the Chebyshev polynomial of degree n. The normal 
computational coordinate r| is related to z through the algebraic 
transformation 

1 + TJ 


z = a- , 

i 2a 

1 + T| 

Z °° ( 2 . 6 ) 

where a is a scaling constant used for proper distribution of points within the 
boundary layer and z^ is the location where free stream asymptotic boundary 


conditions (Malik, Zang, Hussaini 1985) are imposed. 

The spatial discretization employs spectral collocation. The collocation 
points for the periodic directions are 


x j = J L x /Kx, j = 0, 1, ..., Kx - 1 (2.7) 

y m = m Ly/K y , m = 0, 1, ..., Ky - 1. (2.8) 

A staggered grid is employed in the normal direction. Velocities are defined 
at the points 
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tl q = cos 


and the pressures at 


% q 

N r 

q = 0,1 N 

(2.9) 

i ( 7i (q + 1/2 )/N ) , 

q = 0,1 N-l. 

(2.10) 


No artificial pressure boundary conditions are therefore needed. The 
momentum equations are imposed at the points given by Eq. (2.9) and the 
continuity at those given by Eq. (2.10). 


In the spectral collocation method, spatial derivatives of u are obtained 
by differentiating the series expansion coefficients u (t) determined by 
discrete Fourier and Chebyshev transform of the grid-point values of u. The 
temporal discretization involves Crank-Nicolson on the pressure gradient 
and vertical diffusion terms. The remaining terms in the momentum 
equations are handled explicitly by the second-order-accurate Adams- 
Bashforth method. The incompressibility constraint is imposed implicitly. 
The resulting implicit equations are solved iteratively by a minimum 
residual method. The details of the procedure can be found in Malik, Zang 
and Hussaini (1985) and Canuto, Hussaini, Quarteroni and Zang (1988). 

Flow over a curved wall 

The numerical procedure outlined above is applicable to the study of a 
transition problem in two- and three-dimensional boundary layers where 
TS, Gortler or crossflow instability mechanisms may be operative. Here we 
describe the application of the method to the flow over a concavely curved 
plate where both Gortler and TS waves might exist. Extension to three- 
dimensional boundary layers is straightforward. 

Consider an incompressible flow with free stream velocity U,*, along a 

mildly curved wall with constant curvature k = i/r where £ is a characteristic 
length scale and r is the radius of curvature of the wall. If x is the distance 
along the curved wall, y along the span, and z normal to the wall, then the 
governing equations (2. 1-2. 2) may be represented in body-oriented coordinates 
as 

du , du du du - dp 

— + Au — + v — + w — + yuw = - A — 

dt dx dy dz dx 


l 
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( 2 . 11 ) 


2 du 2 . dw 

+ v I V u + y y u + 2yX — 

dz dx 


dv 

+ Xu 

dv 

dv 

dv 

dp 
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dv 
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dx 
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dz 

dw 
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dw 
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dp 




+ V 

+ w 

- y u = - 



dt 


dx 

dy 

dz 


dz 



( 2 . 12 ) 


+ v 


2 dw 2 

V w + y - Y w - 2yX 

dz 


du 

dx 


(2.13) 


where 


du dv 
— + — + 

dx dy 


dw A 

+ yw = 0. 

dz 


(2.14) 


2 2 
2 d d 
V = X — - + — - 

dx dy 
X = ~ ~ — , Y = 


dz 

k 


1 + KZ ’ 1 + KZ 

In equations (2.11-2.14), all velocities are scaled by Uoo , lengths by > l, 


time by MJoo , pressure by pU^; v = l/R, where R is the Reynolds number. A 


Gortler number can be defined as G = rV I k I . 

According to first-order boundary layer theory, the above equations 
reduce to the Blasius flow which provides the basic flow for our study of 
Gortler/TS interaction and transition simulation in a two-dimensional 
boundary layer. 


Linear Stability Theory 

We consider an infinitesimally small disturbance superimposed on 
the Blasius flow. For instance, the streamwise velocity u(x,y,z,t) may be 
written as 
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, j\ tt / \ */ \ i(«x + (3y - cot) 

u(x,y,z,t) =Ug(z) + eu(z)e (2 15) 

under the quasi-parallel flow assumption. Here e is a small parameter, a 
and P are disturbance wave numbers in the x,y directions respectively and 
© = © r + h»i is complex frequency. The real part co r is the actual 
disturbance frequency and the imaginary part coi is the temporal growth 
rate. If coi > 0 the disturbance amplitude increases; otherwise, the 
disturbance decays. 

Linear stability equations may be derived by substituting u(x,y,z,t) 
from equation (2.15) and similar expressions for other dependent variables 
into the Navier-Stokes equations (2.11-2.14) and then retaining only terms of 
0(e). The resulting sixth-order system of ordinary differential equations 
describes an eigenvalue problem for parameters a, P and co which is solved 
by a fourth-order accurate compact difference scheme of Malik, Chuang 
and Hussaini (1982). Depending upon the values of the Gortler number G 
and wave numbers a, P, the governing equations yield a solution for both 
Gortler and TS waves. These solutions provide the initial conditions for the 
nonlinear calculations. In figures (1-2), described below, we provide linear 
results for Gortler and TS waves to emphasize the differences between the 
two instabilities. 

Nonlinear Simulations 

Various nonlinear problems may now be formulated. Suppose, we 
are interested in studying the nonlinear development of Gortler vortices. 
We solve the linear stability problem for the Blasius boundary layer and find 
P = Pg, the growth rate coi and the disturbance eigenfunctions uq, vq, w g 

These are presented in figure 1. The plots show that the stream wise 

/s 

component u g is an order of magnitude larger than the cross-stream 
components v g and w g The spanwise component is 90° out of phase with 

/s /\ 

ug and vq This allows the steady Gortler vortex problem to be represented 

as a system of real equations. In the present case, our complex system of 
equations automatically yields the real solution for the steady Gortler 

/N ^ 

vortex. The structure of v g and Wq is such that it results in a pair of 
counter-rotating vortices . There are higher eigenstates for the Gortler 
problem which yield multiple sets of counter-rotating vortices stacked on 
top of each other. How these higher modes would interact with TS mode 
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may be of interest. The linear stability theory also yields unsteady Gortler 
vortices, whose interaction with TS waves may be of significance, 
particularly in the later stages of transition. Here, however, we only 
consider the steady Gortler vortices. Gortler vortices found in experiments 
are usually steady and of first mode type. In the present study, the steady 
Gortler vortex is assigned an initial finite-amplitude e G . Then an initial 


condition, say for the x component of velocity, is of the form 
u(x,y,z,0) = U^z) + e G Real u G (z)e lpGy ^ 


(2.16) 


Equations (2.11-2.14) are solved with these initial conditions and the 
solution is marched in time. The solution yields the development of the 
fundamental and various harmonics depending upon the number of modes 
allowed in the simulation according to equation (2.5). In the linear regime, 
the temporal development is related to the spatial development through a 
group velocity transformation (c g =dco/dctJ as for TS waves. 

The eigenfunctions for a TS wave are presented in figure 2. The 
eigenfunctions are now complex with disturbance peak lying closer to the 
solid boundary. The interaction of a Gortler vortex with two oblique TS 
waves may be studied by choosing the following initial conditions: 


u(x,y,z,0) = Urfz) + e G Real 


uc(z)e 


iPcy 



«TSX + PTSy 


+ Real 


u(z) e 


i(a<rsx-PTsyj | 


(2.17) 


Various values of PTS in relation to PG may be assigned. For 
example, PTS = PG/2 describes a subharmonic resonance and pTS = PGS 
describes a fundamental resonance. Relative magnitudes of £q and e^g 

describe various regimes of interaction where either the Gortler vortex or 
TS wave is small or when both have large amplitude and interact 
nonlinearly. 

In the Navier-Stokes solution, the mean streamwise velocity U is 
allowed to develop in time according to 
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(2.18) 


au a 2 u _ 

— = — 2 + F > 

at 3z 

U(x,y,z,0) = UB(z). 

Three different forms for the forcing term F were tried: 


2 

„ 3U b 

F = -v — 

dz 


(2.19) 

F = 0, 


(2.20) 

2 

„ dU B 

au B 



> 

(2.21a) 

dz 

dt 


au B „ au B 

at g ax 

> 



(2.21b) 


where Cg is the disturbance group velocity defined above. 

With (2.19), the mean velocity and the boundary layer thickness does 
not change with time. The second condition (F = 0) allows the boundary layer 
thickness to increase and the long-time solution of (2.18) is given by the error 
function. Wray and Hussaini (1984) used this condition for their transition 
simulation in flat plate boundary layers. The condition (2.21a) also allows the 
boundary layer thickness to increase and is designed to give growth in time 
similar to the growth in space according to (2.21b). Since the solutions 
presented below are for short time only, all these formulations essentially 
gave the same results for the development of disturbances of different 
wavenumbers. Therefore, only the results obtained by using equation (2.19) 
are presented below. 

Use of numerical simulation to study the intermodal interaction may 
be described best by providing the example of subharmonic instability in a 
flat plate boundary layer. This instability which is the result of parametric 
resonance between a two-dimensional TS wave and oblique TS waves with 
streamwise wavelength twice that of the fundamental wave was first 
theoretically studied by Herbert (1984) using Floquet analysis. Here we 
consider the fundamental wave to have a = .2033 and prescribe R = 606. An 
amplitude of .01 (i.e., 1%) is assigned to the two-dimensional TS wave. 
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According to Herbert, a strong subharmonic instability of wave 
numbers ((ot/2), (3) should develop for a wide range of P,(see Fig. 9 in Herbert 
1984). The present calculation is performed by including in the initial 
conditions two oblique modes ( a/2, (3 =± .16968 ). Results of the calculation 
are presented in Fig. 3. In this figure, the natural logarithm of the energy 
history of the modes (a,0) and (a/2, (3) are plotted as a function of time. 
Calculations are carried out to about 3.5 linear time periods of the finite- 
amplitude two-dimensional wave. The solid curves are the linear theory 
results for the corresponding modes. The finite-amplitude two-dimensional 
wave follows its linear growth (coi = .7E-5) curve within the computational 
domain. The subharmonic (a/2, |3) follows the linear theory result for the 
primary subharmonic for some time and then diverges from it indicating 
strong secondary instability. The growth rate at the onset of this instability is 
C 0 j * .0098 which drops to (Oj » .0071 towards the end of the computation due 


to nonlinear effects. This is in good agreement with the prediction of Herbert 
for coi = .00824. The subharmonic secondary instability has been linked to the 
squire mode by Herbert (1984). However, in order to capture the instability 
using numerical simulations, we need not provide the actual eigenfunctions 
for the squire mode in the initial conditions. In our calculations we have 
used the eigenfunctions for the primary subharmonic. In another 
calculation the initial conditions for the (a/2,±p) modes were arbitrarily taken 


to be 

u = 0, 


^ 2 -2(z-l) 

v =z e 


and w 


i_dv 

P 5y 


( 2 . 22 ) 


The results are presented in Figure 4. Again the subharmonic secondary 
instability emerges. Most interestingly, both calculations converge towards 
the same growth rate for the subharmonic. Higher harmonics also emerge 
in these nonlinear calculations. In a transition simulation, one would like 
to assign a small random distribution of energy to all the modes, except, of 
course, the finite-amplitude disturbance and let the dominant pattern 
evolve by itself. Such calculations for a flat plate boundary layer have been 
done by Spalart (1984). 
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3. RESULTS FOR GORTLER/TS INTERACTION 

Subharmonic Instability 


Using the multiple scale method, Nayfeh (1981) and Srivastava and 
Dallmann (1987) presented results for Gortler/TS interaction in a boundary 
layer. Their results (derived from Table 3 of Nayfeh (1981) and Figure 5 of 
Srivastava and Dallmann 1987) for R = 950, P = .3 and G = 14 are given in 
Figure 5 where the growth rate of the excited TS waves (in the presence of 
Gortler vortex with 1 % amplitude) of various frequencies is plotted. The 
growth rate of the excited TS waves are shown to increase dramatically in the 
presence of Gortler vortices. It was shown in Srivastava and Dallmann 
(1987) that the growth rate of the excited wave increases with decreasing 
Gortler number and increasing wave number (3. The massive destabilization 
of the TS waves in the presence of Gortler vortices (with only 1 percent 
amplitude) shown in Figure 5 is a bit surprising and of enormous concern to 
a designer. Here we choose this case (R = 950, G = 14, p = .3) for a direct 

numerical simulation. 

The computed growth rate coi for the aforesaid case is .0033 which 
corresponds to a spatial growth rate of a = Rcoi/Cg of 4.48 (Cg was computed 
for the Gortler vortex to be about 0.7). This is in reasonable agreement with 
the spatial calculations of Ragab and Nayfeh (1981). In the nonlinear 
simulations, the Gortler vortex is assigned an initial amplitude of 1%. The 
streamwise wavenumber a was taken to be .103. The complex frequency for 
the oblique (a,p/2) mode is (.038, .00079) and we assign these oblique waves an 
initial amplitude of .1%. The calculation is performed by using 8 Fourier 
modes in the streamwise direction, 16 modes in the spanwise direction and 
33 Chebyshev polynomials across the boundary layer. Occasional tests of 
accuracy were made with twice the number of Fourier modes in streamwise 
and spanwise directions. The time evolution of the energy of six of the modes 
is presented in Fig. 6. The Gortler mode (0,P) continues to grow according to 
its linear theory growth rate. No strong instability in the (a, p/2) mode 
develops at least within the computational period. The growth of this mode 
actually slows down slightly towards the end of the computational period. 
Similar growth behavior of the (a,p/2) mode was observed for other values of P 
at eq = .01 and up to eg = .05. The mode that grows the fastest is (0,2p), i.e. the 
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one with half the wavelength of the finite-amplitude Gortler vortex. This 
seems to agree with the observation of Aihara and Koyama (1981) that the 
"... characteristic spanwise wavelength of importance is shown to be a half 
wavelength of Gortler vortices in the present experiment." Apparently, it is 
the nonlinear distortion of the Gortler vortex that is observed in Aihara and 
Koyama (1981). This nonlinear distortion of the finite-amplitude Gortler 
vortex is neglected in Nayfeh's theory. In our calculation, we also tried to 
"simulate" Floquet theory by zeroing out all higher harmonics after each 
time step. Even this did not result in any change in the growth of the oblique 
TS waves. A careful look at the growth rates given in Nayfeh (1981) for the 
linear Gortler vortex and a comparison with Ragab and Nayfeh (1981) 
indicated that these growth rates could have been obtained only when the 

length scale £ = 1.72^/ v oo x ^oo is used. The length scale used in Nayfeh 

(1981) for T-S problem was £ = V v^x/U^ it was noted by Malik (1986) that 

this discrepancy in the length scales may have been the cause for the 
erroneous growth rates computed by Nayfeh (1981). Indeed, this turned out 
to be the case. Nayfeh and Al-Maaitah (1987) later corrected the length scales 
and found that the multiple scale analysis does not yield any excitation of the 
TS waves when the amplitude of the Gortler vortex is 1 percent. It is still not 
clear why the results of Srivastava and Dallmann (1987) agree with Nayfeh 
(1981). Srivastava and Dallmann also considered combination resonance and 
showed that it is a much more powerful instability mechanism than the 
subharmonic parametric resonance predicted by Nayfeh. They showed that 
the former mechanism produces growth rates which are an order of 
magnitude higher than the growth rates produced by the latter mechanism 
given in Figure 5. However, in view of the disagreement between their 
results and the present simulations of the subharmonic instability, the 
results of Srivastava and Dallmann for the combination resonance should be 
questioned. In any case, results presented by Nayfeh (1981) in Table 1-3 and 
by Srivastava and Dallmann (1987) in Figures 1, 3 and 5 are erroneous. 

The results for Eq = .1 for the above test case are presented in Figure 7. 

Again, we note that the (a, p/2) mode does not grow faster than its linear 
theory results. So the resonance mechanism suggested by Nayfeh (1981) is 
not operative, at least, at this wave number and amplitude of the Gortler 
vortex. Similar calculations using (a,±P) modes in the initial conditions did 
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not indicate any resonant fundamental secondary instability. Later, we will 
show that the (a,(3) mode is excited in the presence of a two-dimensional wave 
((a,0) mode). Another calculation for the subharmonic instability is 
performed at this amplitude (eq = .1) using G = 12 and Pq = .15. The 
stream wise wavenumber for the oblique TS wave is now a = .18 and R = 950. 
The value of co for this TS wave (a, p/2) is co = (.0628, -.0002). So, this wave 
decays linearly. The results of the numerical simulation are presented in 
Figure 8. Now the oblique TS wave is unstable in qualitative agreement with 
the results of Nayfeh and Al-Maaitah (1987) using Floquet theory and 
multiple scale analysis. However, the growth rates are no more than the 
growth rates of the primary Gortler instability. Apparently, the excitation of 
the subharmonic is more pronounced at smaller wavenumbers of the Gortler 
vortex for a given amplitude. This is also in direct contradiction with the 
results of Srivastava and Dallmann who showed that the growth rates of the 
excited oblique wave increase with wave number. In any case, for some 
values of the parameters P and £q, the subharmonic resonance may be 

realized in the Gortler vortex breakdown. Evidently, the Gortler vortex 
amplitude needed for this to happen is in excess of 10%. 


Other Tvnes of Secondary Instabilities 

Nonlinear development of Gortler vortices results in inflectional 
streamwise velocity profiles in the upwash region (peak plane) of the Gortler 
vortex flow. The inflection points caused due to the presence of low velocity 
fluid in the near wall region of the peak planes appear both in the normal 
distribution u(z) and spanwise distribution u(y) (see, e.g., Swearingen and 
Blackwelder 1987). Localized inviscid secondary instabilities may be caused 
due to these inflectional profiles. In the present instance, we use our 
numerical model to investigate possible excitation of two-dimensional (a,0) 
disturbances in the flow field modulated by the presence of Gortler vortices. 
Calculation is first made for the same Gortler vortex as in the example of 
Figure 8. The TS wave now has a = .18 and the linear-theory eigenvalue co = 
(.0613, .0009). The results are presented in Figure 9. The presence of the 
Gortler vortex causes the growth rate of the two-dimensional TS wave to 
increase above its linear theory value. The growth rate of the excited two- 
dimensional TS wave is comparable with the oblique subharmonic instability 
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discussed above. Also shown in the figure are the computed growth history 
of the (a,(5) mode. As noted before, this mode does not get excited in the 
absence of the (a,0) mode. However, now the (a,P) mode attains higher 
growth rates in the presence of the (a,0) mode. The growth rate of this mode 
is compared with other modes in Figure 10 at time t - 400 (which is 
equivalent to about 4 periods of the two-dimensional mode) along with the 
corresponding linear theory results. The figure also contains growth rate 
results for (a, (3/4) mode which is also excited. 

Another calculation is performed at G = 14, (3q = .5 with the same 
amplitude for the Gortler vortex as before (eq. = .1). The two-dimensional TS 
wave now has a = .15 and co = (.04940, .00227). The results are presented in 
Figure 11. Again, the two-dimensional TS wave attains a growth-rate 
much higher than its linear value. The (a,P) mode also grows fast due to 
nonlinear interaction. The orientation \|/ = tan* 1 (3/a for this mode is 73.3°. 
The orientation of the (a,P) mode of Figure 11 was 39.8°. So there seems to 
be a wide range of angles for which the (a,P) mode is unstable. In this 
example, the two-dimensional TS wave was linearly unstable. Another 
calculation is presented in Figure 12 where the two-dimensional TS wave (a 
= .2) is linearly stable (linear theory co = (.0693, -.00115)). As shown in the 
figure, both the (a,0) and (a,P) modes grow in the nonlinear calculations 
regardless of the fact that the TS mode chosen is linearly stable. This is 
consistent with the results presented in Figures 3 and 4 where it was 
shown that the form of the initial conditions for the excited mode do not play 
a critical role. As noted previously, the first harmonic of the Gortler vortex 
gains energy faster than any other mode indicating a spectral broadening 
process. However, the three-dimensionality and unsteadiness is caused by 
the development of (a,0), (a,P) or possibly (a, p/2) modes depending upon the 
initial conditions and flow parameters involved. The computed growth rate 
of these excited modes, at least initially, remain small when compared to 
the growth rate of the Gortler vortex itself. However, the Gortler vortex does 
reach a saturated state (as also noted in the experiments of Swearingen and 
Blackwelder 1987) at which time the secondary instability modes will be 
expected to grow much faster and would lead to various kinds of breakdown 
mechanisms as observed by Swearingen and Blackwelder (1987), Bippes 
(1978) and Aihara and Koyama (1981). Calculations indicate that at 
saturation, the Gortler vortex amplitudes are in excess of 20 percent. We 
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believe that our numerical model would be able to provide detailed 
qualitative information about the development of unsteadiness and eventual 
breakdown; however, this calculation would require much higher 
resolution than used here, especially in the normal direction in order to 
resolve critical layers that are present near the edge of the boundary layer. 

4. CONCLUDING REMARKS 

Direct numerical simulation of wave-interactions is useful not only 
in testing theories but also in bringing out new instability mechanisms in 
boundary layers. In this paper, we have presented a numerical model for 
Gortler/TS interaction which may be used to investigate various possible 
regimes of interactions between Gortler vortices and Tollmien-Schlichting 
waves. In the numerical examples, we have only considered the regime 
where Gortler vortices have finite-amplitude and TS waves are relatively 
small. Even though our model is based upon the parallel flow 
approximation, it provides the qualitative features of the possible 
interaction mechanisms. Of course, any mechanism that is mainly caused 
by non-parallel effects cannot be treated by our model. What mechanisms 
might be ruled out because of this limitation is open for discussion. 
Gortler/TS interaction takes place at relatively large amplitudes of Gortler 
vortices (as primary instability) and therefore secondary instability 
mechanisms leading to unsteadiness most likely involve nonlinearities, 
and the direct numerical simulation can obviously take account of such 
nonlinearities. Finally, we have found that the Gortler/TS interaction 
results presented by Nayfeh (1981) and Srivastava and Dallmann (1987) are 
erroneous. 

Acknowledgements 

This work was partly supported by NASA Contracts NAS 1-18240 
(MRM) and NAS1-18107 (MYH). 


14 


REFERENCES 


Aihara, Y. & Koyama, H. 1981 Trans. Japan Soc. Aero. Space Sci. 24, 78. 
Bippes, H. 1978 NASA TM- 75243. 

Canuto, C., Hussaini, M. Y., Quarteroni, A., & Zang, T. A. 1988 Spectral 
Methods in Fluid Dynamics, Springer-Verlag, New York. 

Floryan, J. M. & Saric, W. S. 1980 AIAA Paper No. 80-1376. 

Hall, P. & Smith, F. T. 1988 Proc. Roy. Soc. Lond. A417, 244. 

Herbert, T. 1984 AIAA Paper No. 84-0009. 

Herbert, T. & Morkovin, M. V. 1979 Laminar-Turbulent Transition, 37. 
Malik, M. R. 1986 AIAA Paper No. 86-1129. 

Malik, M. R., Chuang, S., & Hussaini, M. Y. 1982 ZAMP 33, 189. 

Malik, M. R., Zang, T. A., & Hussaini, M. Y. 1985 J. Comp. Phys. 61, 64. 
Nayfeh, A. H. 1981 J. Fluid Mech. 107, 441. 

Nayfeh, A. H. & Al-Maaitah, A. 1987 AIAA Paper No. 87-1206. 

Ragab, S. A. & Nayfeh, A. H. 1981 The Physics of Fluids 24, 1405. 

Spalart, P. R. 1984 9th Inti. Conf. on Numer. Meth. in Fluid Dynamics. 
Srivastava, K. M. & Dallmann, U. 1987 Phys. Fluids 30, 1005. 

Swearingen, J. D. & Blackwelder, R. F. 1987 J. Fluid Mech., 182, 255. 

Wray, A. & Hussaini, M. Y. 1984 Proc. Roy. Soc. London A 392, 373. 

Zang, T. A. & Hussaini, M. Y. 1985 AIAA Paper No. 85-0296. 


15 


1.2 



Figure 1. Eigenfunctions for a Gortler vortex in a blasius boundary layer: 
G = 14, P = .3. 
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layer: R = 950, a = .103, P = .15. 
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Fig. 2(b) imaginary part 
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Figure 3. Computed evolution of energies of a two-dimensional TS wave 
((a,0) mode) and the oblique subharmonic ((ot/2,(3) mode). In this 

calculation, R = 606, a = .2033, P = .16968. A strong secondary 
subharmonic instability develops. Included in the initial 
conditions are two oblique primary subharmonics. The solid 
lines represent linear theory results. 


19 



12 


ax 10 



FxlO 6 


Figure 5. Spatial growth rates of excited oblique TS waves in the presence 
of a Gortler vortex with 1% amplitude (from Nayfeh 1981 and 
Srivastava and Dallmann 1987) and linear theory results. As we 
show later, the curve labelled "excited TS wave" is erroneous. 
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Figure 6. Computed evolution of energies of various modes in the presence 
of a Gortler vortex ((0,(3) mode) with 1% amplitude in a two- 
dimensional boundary layer on a concave wall. Two oblique TS 
waves ((a,±(3 /2) mode) with amplitudes of .1% are included in 
the initial conditions. In this calculations R = 950, G = 14, a = 
.103, (3 = .3. The solid lines represent linear theory results. 
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Figure 9. The effect of a GQrtler vortex ((0,(3) mode) with 10% initial 
amplitude on a two-dimensional TS wave ((a,0) mode). In this 
calculation, R = 950, G = 12, a = .18, (3 = .15. The (a, (3) mode 
grows due to nonlinear interaction. The solid line represents 
linear theory result. 
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Figure 10. Temporal growth rates of various modes at time t = 400 (a = .18, 
p = .15). 
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